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FOREWORD 


This  report  was  prepared  in  the  Mechanical  Engineering 
Department  of  The  Johns  Hopkins  University  on  Contract  AF- 
33(6l6)-152,  under  Project  No.  3066,  "Gas  Turbine  Technology," 
Task  No.  70153,  "Compressor  Blading  Systems  of  High-Deflection 
Angle.  "  The  report  covers  the  last  phase  of  the  contract  work, 
which  was  under  the  direction  of  Dr.  H.  Yeh.  Previous  phases  of 
this  contract  work,  which  were  under  the  direction  of  Dr.  G.  Wis- 
licenus,  resulted  in  the  preparation  of  WADC  TR  55-348,  "Three- 
Dimensional  Flow  in  Turbomachinery,  "  Volumes  I  and  II,  by 
Leroy  Smith. 

This  contract  work  was  administered  by  the  Fluid  Dynamics 
Research  Branch  of  the  Aeronautical  Research  Laboratory, 

Wright  Air  Development  Center,  with  Mr.  S.  Hasinger  acting  as 
Task  Scientist. 
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ABSTRACT 


Fluctuation  intensities  at  selected  stations  in  am  axial-flow 
compressor  stage  were  measured.  These  measurements  recorded 
the  total  fluctuation  which  consists  of  both  periodic  amd  random 
(turbulence)  parts.  Through  the  investigation  on  the  periodic  part 
of  the  fluctuation,  an  estimate  can  be  made  on  the  turbulence  in¬ 
tensity. 

The  intensity  of  the  induced  velocity  due  to  the  rotor  at  points 
in  front  of  the  rotor  is  analyzed.  Also  investigated  is  the  fluctuation 
intensity  due  to  the  periodic  presence  of  wakes  behind  the  rotor. 

Such  analytical  considerations  are  found  to  be  in  substantial  agree¬ 
ment  with  the  measured  data. 

The  wake  behind  blades  is  found  to  decay  and  spread  in  the 
same  manner  as  the  wake  behind  single  airfoils  except  in  the  im¬ 
mediate  neighborhood  of  the  following  blade  row,  at  which  place  the 
decay  is  accelerated. 
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1.  Introduction 


Recent  interest  on  the  investigation  of  nonsteady  flows  in  an 
axial-flow  turbomachine  seems  to  be  centered  mostly  on  the  phenom¬ 
enon  of  stall  propagation  and  surge.  Results  on  this  type  of  measure¬ 
ment  have  been  reported  by  Huppert  and  Benser  (Ref.  1),  Iura  and 
Rannie  (Ref.  2),  and  Emmons,  Pearson  and  Grant  (Ref.  3).  By  com¬ 
parison,  few  investigations  on  the  nonsteady  flows  in  the  normal 
operating  range  of  an  axial-flow  turbomachine  have  been  reported. 

This  is  undoubtedly  due  to  the  general  belief  that  as  long  as  the 
machine  is  not  in  the  stalling  range,  the  nonsteady  part  of  the  flow 
due  to  turbulence  and/or  the  relative  motion  between  neighboring  blade 
rows  does  not  influence  the  performance  of  the  machine  to  any  great 
extent.  Hence  present  day  design  is  still  based  on  the  assumption  of 
steadiness  and  often  also  of  axial- symm etry,  modified  and  refined  if 
necessary  by  the  consideration  of  (axially  non-symmetrical)  secondary 
flows.  (Secondary  flow  in  its  commonly  accepted  sense  is  still  a  steady 
flow  phenomenon  relative  to  the  blade  passage.  )  A  previous  report  by 
Smith  (Ref.  4)  issued  in  this  laboratory  under  the  same  contract  is  one 
example  of  the  effort  along  that  direction. 

It  is  inevitable,  however,  that  as  the  art  of  the  science  ad¬ 
vances,  a  deeper  insight  into  the  complicated  flow  phenomenon  in  an 

NOTE:  This  report  was  released  for  publication  in  December  1955. 
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axial-flow  turbomachine  can  be  gained  only  through  the  additional 
investigation  of  the  non- steady  part  of  the  flow.  Kemp  and  Sears 
(Ref.  5)  were  the  first  ones  to  report  on  the  influence  on  the  overall 
characteristics  (such  as  lift,  moment,  and  induced  drag)  of  the  blades 
of  a  cascade  due  to  the  interaction  between  stationary  and  moving  blade 
rows.  A  different  formulation  on  this  problem  by  means  of  conformal 
transformation  was  reported  by  Woods  (Ref.  6).  Recently  in  this 
laboratory,  Meyer  (Ref.  7)  made  a  detailed  calculation  on  the  nonsteady 
parts  of  pressure  and  velocity  distribution  on  a  blade  surface  due  to 
the  action  of  the  travelling  viscous  wakes  shed  from  an  upstream  cas¬ 
cade  of  blades.  This  type  of  flow,  although  nonsteady,  is  still  periodic. 
On  the  other  hand,  there  exists  also  the  random  type  of  nonsteady  flow 
or  turbulence.  It  is  to  be  expected  that  either  type  of  fluctuations  will 
effect  the  profile  drag  (as  distinct  from  the  induced  drag)  and  the 
maximum  attainable  lift  coefficient  of  the  blades.  Furthermore,  in 
at  least  one  case  the  non- steady  part  of  the  flow  may  be  decisive  in 
the  normal  performance  of  the  machine.  This  is  the  case  of  com¬ 
pressors  with  small  blades  operating  at  high  altitudes  so  that  Reynolds 
number  is  low  enough  to  be  in  the  transition  range.  Whether  the 
boundary  layer  along  the  blades  will  be  mostly  laminar  or  mostly  tur¬ 
bulent  or  whether  there  will  be  laminar  separation  will  presumably 
depend  on  the  nonsteady  as  well  as  the  steady  character  of  the  flow. 


This  report  presents  a  series  of  preliminary  nonsteady  velocity 
measurements  in  a  single-stage  experimental  compressor  with  Reynolds 
number  approximately  20,  000  (based  on  blade  chord  length).  Results 
include  the  intensity  and  the  decay  of  wakes,  the  induced  velocity  up¬ 
stream  of  the  rotor  and  the  turbulence  level  throughout  the  machine. 
Analysis  was  made  insofar  as  possible  to  interpret  the  experimental 
results. 


2-a.  Compressor  (Fig.  1).  The  low-speed  single-stage  com¬ 
pressor  used  in  these  experiments  has  been  fully  described  in  a  pre¬ 
vious  report  (Ref.  4),  Briefly,  the  unit  consists  of  a  row  of  inlet 
guide  vanes,  followed  by  a  rotor  and  a  stator.  It  was  designed  so  that 
the  energy  addition  is  proportional  to  the  radius.  The  hub-tip  diameter 
ratio  is  0.5.  A  cross-sectional  view  of  the  unit  indicating  all  the 
major  dimensions  is  shown  in  Fig.  1. 

One  alteration  on  the  unit  has  been  made  since  the  previous 
report.  In  order  to  investigate  fully  the  decay  and'  the  spread  of  the 
wake  behind  the  inlet  guide  vanes,  the  distance  between  the  inlet  guide 
vanes  and  the  rotor  has  been  increased  by  the  addition  of  a  spacer 
section  of  5-1/16"  in  axial  length.  Fig.  1  shows  this  altered 
dimension. 
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2-b.  Traversing  Apparatus  (Fig.  2).  The  holder  for  the  con¬ 
ventional  (static  pressure,  total  pressure,  3-claw  yaw,  etc.  )  as  well 
as  the  hot-wire  probes  is  supported  by  two  perpendicular  sets  of  rails. 
The  whole  unit  is  designed  in  such  a  way  that  complete  traversing  of 
the  probe  can  be  made  along  axial,  circumferential  and  radial  direc¬ 
tions.  The  probe  can  also  be  rotated  about  its  own  axis. 

In  order  to  make  a  detailed  survey  of  the  wake  behind  the 
inlet  guide  vanes,  it,  is  necessary  to  be  able  to  insert  the  probe  at  a 
large  number  of  stations  behind  the  vanes  in  both  tangential  and  axial 
directions.  To  avoid  the  drilling  of  a  large  number  of  holes  in  the 
casing,  a  portion  (4"  X  10")  of  the  casing  wall  is  replaced  by  a  lattice 
of  interchangeable  sliding  plexiglas  blocks  (Fig.  2).  One  of  the  blocks 
is  provided  with  holes  at  desired  locations  through  which  a  probe  can 
be  inserted.  This  "master"  block  can  be  shifted  to  any  tangential  and 
axial  station  and  the  entire  lattice  remains  airtight.  A  circumferential 
survey  at  intervals  of  1  degree  can  be  easily  made.  This  device  proves 
to  be  so  convenient  that  the  authors  would  like  to  recommend  it  for 
future  such  surveys. 

2-c.  Hot-wire  Anemometer  (Fig.  3).  To  measure  velocity 
fluctuations,  a  tungsten  wire  of  .  0002"  diameter  and  about  1/16"  long 
was  used.  The  cold  resistance  of  the  wire  was  approximately  6  ohms, 
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and  its  time  constant  about  0.  5  /U  sec.  The  ratio  of  hot  to  cold  re¬ 
sistance  was  fixed  at  1.  45. 

The  hot-wire  anemometer  equipment  used  in  these  experiments 
is  of  the  constant-current  type  based  on  the  design  by  Kovasznay 
(Ref.  8).  The  first  stage  of  the  amplifier  has  been  modified  to  further 
reduce  the  noise  level.  A  resistance-capacitance  circuit  is  used  for 
full  compensation  of  the  time  lag  of  the  wire.  The  compensating 
amplifier  gives  essentially  flat  response  up  to  10,  000  cps.  The  noise 
level  at  the  operating  condition  of  these  experiments  is  less  than  0.1  mv 
at  input.  This  corresponds  to  a  turbulence  intensity  of  much  less  than 
0.  001  (0. 1%)  in  these  measurements. 


3.  Description  of  Experiments 


All  measurements  described  in  this  report  were  made  at  one 
operating  condition  which  was  near  the  design  point  of  the  compressor 
and  well  to  the  right  of  the  peak  in  the  head-capacity  curve.  Hence 
the  blades  were  not  near  stalling.  For  these  preliminary  measure¬ 
ments,  the  radial  location  of  all  probe  points  was  arbitrarily  fixed  at 
4-3/16"  from  the  blade  tip.  This  was  roughly  near  the  middle  span 
of  the  blade. 

The  measuring  stations  along  the  axial  direction  were  as 
follows  (in  the  order  of  direction  of  flow): 
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Before  inlet  guide  vanes  -  Station  1 

Between  inlet  guide  vanes  and  rotor  -  Stations  2,  3,  4,  5,  6 
Between  rotor  and  stator  -  Stations  7,  8 
Between  stator  and  exit  -  Stations  9,  10 
The  axial  distances  of  these  stations  are  shown  in  Fig.  1  and  repeated 
in  Fig.  4.  The  latter  figure  is  a  cylindrical  sectional  view  at  the 
measuring  radius. 

•  4.  Results  of  Experiments 

A  detailed  survey  of  the  flow  behind  the  inlet  guide  vanes  was 
made  at  Stations  2  to  6.  The  measurements  included  direction  of 
flow  (by  three-claw  yaw  probe),  total  pressure  (by  Kiel  probe)  and 
static  pressure  (by  spherical  probe).  Distributions  of  dynamic  pressure 
at  various  s.tationsof  the  wake  behind  inlet  guide  vanes  are  shown  in 
Fig.  5. 

The  intensity  of  velocity  fluctuation  n/u^  /U  (defined  as  the 
ratio  of  the  r.  m.  s.  of  the  fluctuating  velocity  u  to  the  mean 
velocity  U)  at  all  stations  is  as  follows.  (Unless  stated  otherwise, 
the  hot-wire  is  normal  to  the  mean  direction  of  flow.  ) 
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Fluctuation  Intensity 

Location 

Station 

Case 

nTu2  /U  (%) 

Before  I.  G.  V. 

1 

1 

0.  24 

r  2 

2a 

in  middle  of  wake 

6.  13 

2b 

outside  of  wake 

.  39 

3 

3a 

in  middle  of  wake 

3.  24 

3b 

outside  of  wake 

.35 

1  4 

4a 

in  middle  of  wake 

.  3.  00 

Between  I.  G.  V.  < 

4b 

outside  of  wake 

.  30 

and  rotor 

5 

5a 

in  middle  of  wake 

2.  84 

5b 

outside  of  wake 

.  45 

^  6 

6a 

in  middle  of  wake 

2.  93 

6b 

outside  of  wake 

1.  45 

(  7 

7a 

wire  in  tangential  dir. 

5.  47 

Between  rotor 

J 

7b 

wire  in  axial  dir. 

2.  20 

and  stator 

1  a 

8 

wire  in  tangential  dir. 

3. 16 

f  •> 

9 

wire  in  tangential  dir. 

1.30 

Behind  stator  ' 

l  .0 

10 

wire  in  tangential  dir. 

1.  32 

Oscilloscope  photographs  of  the  instantaneous  velocity  for  all 
cases  (except  Cases  4a  and  4b)  are  reproduced  in  Fig.  6  to  11.  The 
faint  detached  lines  on  top  of  the  film  strips  are  the  images  of  a 
small  neon  bulb  which  flashes  120  times  a  second.  These  lines,  there¬ 
fore,  fix  the  time  scale  of  the  abscissa.  A  sine  wave  of  known  voltage 
was  fed  into  the  amplifier  (without  compensation)  and  the  image  on 
the  oscilloscope  screen  was  photographed  at  the  beginning  of  each 
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film  strip.  The  voltage  of  the  sine  wave  was  adjusted  in  each  case  to 
give  a  picture  of  approximately  the  same  height  as  that  of  the  velocity 
fluctuations.  All  controls  in  the  amplifier  and  the  oscilloscope  were 
left  untouched  between  the  taking  of  the  sine  wave  and  that  of  the  in¬ 
stantaneous  velocity.  In  this  way  the  scale  of  the  ordinate  of  the 
film  strips  is  fixed  in  terms  of  input  voltage  of  the  amplifier.  On  the 
other  hand,  the  relation  between  the  voltage  and  the  velocity  fluctua¬ 
tions  are  established  by  the  hot-wire.  We  therefore  have  a  convenient 
means  of  estimating  the  magnitude  of  all  instantaneous  velocities. 


5.  Analysis  and  Interpretation 


5a.  Separation  of  Periodic  and  Random  Fluctuations.  Let 
the  instantaneous  velocity  U  be  the  sum  of  an  average  velocity  U 
and  a  (total)  fluctuating  velocity  u.  The  fluctuating  velocity  u  in 
turn  consists  of  periodic  fluctuations  Up  and  random  fluctuations 
ur.  (We  here  use  the  word  "periodic"  to  describe  the  property  of 
repeating  from  one  rotor  blade  to  another.  )  We  have  thus, 

(J  (t)  =  (J  +  Vtp(t  )  +  U.r(t)  (1) 


By  definition 

lAp  +  vlr  =  0 


(2) 
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where  a  bar  on  top  means  time  average.  Eq.  (2)  indicates  that 
Up  =  -Ur  although  Up  and  ur  may  not  be  zero  individually. 
However,  we  can  easily  remove  this  complication  by  letting 


up(t)  =  up  +  up(t) 


ur(t)  =  ur  +  u'r(t) 


(3) 


Now 


u^  =  *  Uj?  =0  by  definition. 


(4) 


Also  Up  +  ur  =  Up1  +  ur  1  by  virtue  of  Eq.  (2).  The  mean 
square  of  the  total  fluctuations  becomes  thus, 


(Uq  +  ur  )2  =  (Uq  +  u£  )2 


'  2  /  2  ^  /  / 

Up  +  ur  +  2  Up  ur 


p  r '  '  p 

Now  since  ur‘  and  Up*  are  uncorrelated  by  definition,  the  last 


term  2up'  ur'  becomes  2  up '  ur'  (see  for  example  Cramer, 
Ref.  9)  which  is  zero  on  account  of  Eq.  (4).  Hence  we  have 


*/u"z  =  V (Up  +  ur  =  Ju£  +J u'f 


(5) 


Eq.  (5)  gives  the  intensity  of  total  fluctuations  in  terms  of  the  in¬ 
tensities  of  the  individual  periodic  and  random  fluctuations.  The 
experiment  measures  the  first  term,  namely  \TW*  .  In  order  to 
make  use  of  Eq.  (5),  we  have  to  know  in  addition  either  uT  or 


ur'^.  It  is  for  this  reason  that  the  evaluations  of  Up'  before 


and  after  the  rotor  are  presented  in  sections  5c  and  5d, 
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5b.  Wake.  Fig.  5  shows  the  familiar  phenomenon  that  the 
wake  behind  a  body  widens  in  the  direction  of  flow  while  at  the  same 
time  decreases  its  intensity.  It  is  interesting  to  compare  our  case 
of  wake  behind  inlet  guide  vanes  with  that  behind  a  symmetrical  two- 
dimensional  body  for  'tohich  the  following  relations  hold.  (See  Durand, 
Ref.  10,  or  Goldstein,  Ref.  11.  ) 


Ud. 


max 


oC 


yD  OC  x‘/jl 


(6) 

(7) 


where  Uj  =  maximum  velocity  defect  in  the  wake, 

max 

yQ  =  half-width  of  the  wake,  i.  e.  ,  the  distance 

between  the  center  of  the  wake  and  its  edge 
where  U  =  0. 

and  x  =  distance  along  direction  of  flow  measured 

from  a  suitable  origin. 

The  maximum  velocity  defect  L^max  in  these  experiments  can  be 

easily  calculated  from  the  dynamic  pressure  distribution.  Fig.  12 
2 

shows  a  plot  of  versus  distance  x  is  rendered 

non-dimensional  as  a  ratio  to  the  maximum  velocity  outside  of  the 
wake  at  Station  2).  According  to  Eq.  (6)  this  plot  should  be  a 
straight  line  with  the  intercept  on  the  abscissa  at  some  distance 
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1 
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upstream  of  the  trailing  edge  of  the  vanes.  Fig.  12  shows  that  this  is 
indeed  true  for  a  distance  of  about  2  chord  lengths  downstream  of  the 
guide  vanes.  The  intercept  of  the  straight  line  on  the  abscissa  is  at  a 
distance  of  approximately  20%  of  chord  length  upstream  of  the  trailing 
edge.  This  agrees  well  with  the  result  on  single  airfoils  by  Silverstein, 
Katzoff  and  Bullivant  (Kef.  12).  However,  after  a  distance  of  two  chord 
lengths  downstream  the  decay  of  the  wake  is  greatly  accelerated  so 
that  at  Station  6  the  value  for  l/Ud^Jax  is  considerably  more  than 
what  it  should  be  from  the  straight  line  relation.  An  obvious  explana¬ 
tion  is  that  the  nonsteady  induced  velocity  from  the  downstream  rotor 
tends  to  increase  the  mixing  of  the  non-uniform  flow  in  the  wake.  The 


r 

verification  of  this  is  the  fact  that  such  an  induced  velocity  is  found 
to  be  indeed  appreciable  at  the  same  station.  (See  Section  5c.  ) 
j  An  attempt  to  plot  wake  width  against  distance  to  check  Eq.  (7) 

runs  into  the  difficulty  that  the  width  is  not  well-defined.  One  could  of 
course  invent  a  representative  width  such  as,  for  example,  the  ratio 

.  of  the  integrated  area  of  the  (velocity  defect)  distribution  curve 

s' 

to  Udmax*  However,  such  a  procedure  does  not  seem  to  be  worth- 

ji. 

j  while  when  the  distribution  of  dynamic  pressure  is  not  at  all  smooth 

t 


as  can  be  seen  in  Fig.  5.  Nevertheless  a  rough  check  was  made,  and, 
according  to  expectation,  it  was  again  found  that  the  rate  of  spread  of 
wake  follows  £q.  (7)  up  to  about  two  chord-lengths  downstream, 
after  that  the  rate  of  spread  is  accelerated. 
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The  fluctuation  intensity  »>/  u^  /U  between  the  inlet  guide  vanes 


and  the  rotor  is  plotted  against  the  axial  distance  in  Fig.  13.  The  upper 
curve  shows  that  the  maximum  turbulence  intensity  in  the  wake  de¬ 
creases  with  distance  in  the  direction  of  flow.  This  result  is  similar 
to  that  reported  by  Townsend  (Ref.  13,  14)  for  the  wake  behind  a  circular 
cylinder.  Since  Townsend's  measurements  started  at  a  distance  of  80 
times  cylinder  diameter  downstream  of  the  cylinder  whereas  the  last 
station  (i.  e.  ,  Station  6)  in  these  measurements  was  only  at  a  distance 
of  24  times  maximum  blade  thickness  from  the  trailing  edge  of  the 
blades,  no  comparison  between  the  two  cases  is  possible. 

The  lower  curve  in  Figure  13  shows  that  the  fluctuation  in¬ 
tensity  outside  of  the  wakes  remains  at  a  constantly  low  level  for  a 
distance  of  approximately  two  chord  lengths  behind  the  inlet  guide 
vanes,  after  that  the  intensity  increases  rapidly  due  to  the  periodic, 
induced  velocity  of  the  downstream  rotor.  Examination  of  oscillo¬ 
scope  photographs  shows  that  at  Station  5  outside  of  wakes  (Fig.  8, 

Case  5b)  the  nonsteady  velocity  begins  to  show  signs  of  periodicity 
due  to  the  passing  of  the  rotor.  At  Station  6  (Fig.  9,  Case  6b)  the 
periodic  velocity  completely  overwhelms  the  random  velocity  or  tur¬ 
bulence.  On  the  other  hand,  the  nonsteady  velocity  in  the  wake  at  the 
same  stations  is  still  predominated  by  turbulence,  since  the  amount 
of  periodic  velocity  is  much  less  than  the  amount  of  turbulence  there. 

The  following  section  is  an  attempt  to  estimate  the  magnitude 
of  the  induced  velocities  such  as  the  one  shown  in  Fig.  9,  Case  6b. 
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5c.  Induced  Velocities  due  to  a  Moving  Cascade.  To  analyze 
the  nonsteady  velocities  induced  upstream  due  to  the  motion  of  a  rotor 


downstream,  it  is  necessary  to  simplify  the  problem  by  considering 
the  rotor  as  a  plane  cascade  of  infinitely  many  blades  or  airfoils 
bounded  by  two  walls  representing  the  casing  and  the  hub  (Fig.  14).  We 


shall,  in  this  section,  use  the  following  notation  (see  Fig.  14): 

x  -  direction  along  the  span  of  the  airfoils,  or  radial 
direction  in  an  axial-flow  turbomachine. 

y  -  direction  of  the  "axis"  of  the  cascade,  or  tangential 
direction  in  an  axial-flow  turbomachine. 

z  -  direction  perpendicular  to  x  and  y,  pr  axial 
direction  in  an  axial-flow  turbomachine. 

^  -  direction  of  the  chord  of  the  airfoils.  The  angle 

between  z  and  §  is  denoted  by  ot  . 

7J  -  direction  perpendicular  to  £  ,  in  the  y-z  plane. 

u,  v,  w  -  induced  velocity  components  along  x,  y,  z 
respectively. 


-  induced  velocity  along  • 

Considering  one  airfoil  of  the  cascade  (Fig.  15),  its  action  can 
be  approximated  by  a  series  of  vortex  lines  (bound  vortices)  whose 


circulation  T  is  in  general  a  variable  along  both  the  chord  and  the 
span  of  the  blade.  Thus 

r  =  y  <  §  ,  x  ) 
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Because  the  circulation  varies  along  the  vortex  line,  trailing 
vortices  will  be  shed.  According  to  Helmholtz'  Law,  the  strength 

^  y 

of  the  trailing  vortices  is  equal  to  ~^x  •  ^*or  ^rst  approxima¬ 

tion,  it  is  assumed  that  the  trailing  vortices  move  downstream 
with  the  main  velocity  which,  for  cascades  of  thin  airfoils,  has 
the  direction  of  $  .  Since  ~Y  varies  continuously  along  x 

for  a  bound  vortex  line  stationed  at  £  =  _§  ,  the  resulting 

trailing  vortex  lines  will  smear  over  the  entire  £  -x  plane, 
from  §  =  3  t°  £  =  °°  •  The  same  holds  true  for  all 
other  bound  vortices.  Finally,  the  blade  is  bounded  by  an  upper 
and  a  lower  wall  (casing  and  hub)  whose  action  is  identical  to  having 
the  same  vortex  pattern  repeated  as  images  to  the  two  walls.  Thus, 
if  [_  is  the  blade  length,  we  must  have 

Yr-x)  =  Y  y  fL  +X)  =  Y(L-X) 

This  is  achieved  by  expressing  Y"  in  a  Fourier  cosine  series, 

v  a  x)  =  x  f\„  (s)  (8) 

n-o 

(1)  Bound  Vortices.  Consider  the  induced  velocity  at  any 
upstream  point  P  (xQ  ,  -yQ  ,  -z0  )  in  the  flow  field  due  to  a  single 

bound  vortex  line  extending  from  x  =  -o  to  x  =  +0  and  located  at 
y  =  y,  z  =  z  (Fig.  16).  The  induced  velocities  along  y- and  z- 


1 1 


+ 

•  4». 


directions  due  to  such  a  single  vortex  line  are: 


A»t> 


V  = 


•oo 


v  (x,  l)  (Za  +  fcAjcx.)  <tx _ 

47T  [(%  +  fSi»«)*+(k,  +f 


„+°° 


ur  St 


Y  (X,  j )(y<,  +  f  S//tet )  c/x 


-  oo 


^(O/o  +fS"»*)2-i  ( Z0+  ?c°soc)2-i-  CX-X^JH 
a 

Integrating  through  the  entire  chord  from  §  =  0  to 
J-  =  c  (see  Fig.  16)  and  summing  up  all  blades  in  the  cascade, 
we  obtain  the  induced  velocities  due  to  all  the  bound  vortices  repre¬ 
senting  the  infinite  cascade: 
oo  .c  j-oo 


V 


-w 

n*°i  J.oo 


Y  (  x,  V  (Zo  i-J  «)  <*f 


oo 


W 


-i 


f‘7r£(%+£s'»*  ± ns)z+  ( £0  +j&>to()2+  (x-Xo)2^^2- 

<’> 

Y(X,  f)(X  +JSAnt  ±  os)  <U.  <*£ 


n*°Jo  J  +fjt»o(± /?s)2+  (Z,  +  J <*s<xf+(x~Xo)*l*$. 

—  OO  — * 


The  evaluation  of  the  above  integrals  in  the  general  case  will  be 
quite  difficult.  It  can  be  seen  from  Eq.  (8)  for  y  (x,  f  )  that 
the  computation  would  be  facilitated  if  one  had  a  table  of  the  follow¬ 
ing  definite  integral: 

+  oo 


I  («.,  Xo)  -  f  £ ffe _ 

a*  +  (X-x.)2JsYi 


Such  a  table,  however,  does  not  seem  to  be  available. 
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We  consider  next  the  following  special  case.  If  Y  -  ~Y  {  ) 

only,  Eqs.  (9)  then  become 

Oo  z'  C 


y  f  ($)(?..+£  «)  d£ _ 

n=o  J0  2  7T[(yo  +  SSin  oi  ±  hs)2+( +S&so()zJ 

J-  Y(Z)  CX  +  fS/no<  4  hs)  df _ 

7  =  0  -o  2n  [  ( y.  +§s/n*  ins')*  +  (Z,-f  fcos&f'j 


oo  C 


/?  =  <>  0 


To  reduce  the  above  equations,  one  makes  use  of  the  fact  that 
v  and  -w  are  the  real  and  the  imaginary  parts  of  the  following 
function, 


]r  fC  -r  (S)  d£ _ 

>=*  Jt  Z7r(i(Z  +  ±ns)  +  (ZB  +  J c*s<x )J 


V  -  L  ur  = 


This  equation  can  obviously  be  obtained  also  by  considering  the 
complex  potential  of  two-dimensional  vortices  to  which  the  problem 
is  now  reduced. 


From  the  relation 


-f  oo 

coM  (x+iy)  «  r  _ L 

7  ~  +  i 


(  y-tyvj 


we  can  write 


T  _ L _ _ 

„=0  ( 2.  «)  +  i  (X  +  S  Sin*  -±  ns) 

=  f  ■=»■*  ffe.+- iy.  +fel«) 
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Hence 


V-  i  us 


2  5 


j1  f*.  *  +f  *L*) y(S)  <*£ 


(ii) 


We  can  now  separate  the  real  and  imaginary  parts,  since 

Coih  (x+  iy)  -  X  -  t  -/'guff  Sec/)  Z 2 

~f~anft  x  ~f~  ~htn2 y 

we  have 

_  j_  f  ta*L + fCcs go Sec*?-(' & thii *L  v/ri  ^ 

2  5  Ja  fCnsct)  +  t ‘an2?! (t/o  -f-fS/nvt) 

c  (12) 
ur«=  ~L  f -^tCX  +?*'">«} Secf>2:r(%0  +  Sc*s«.)  ydjjf 

2  s  Jo  ^2f(£a  +*•*«)+  ■/*» 


For  the  special  case  of  V  =  const,  from  Jj"  =  0  to  J  =  c 
(rectangular  loading),  Eq.  (11)  can  be  integrated  as  follows, 


V-  i  MJ  —  — 
271 


e 


3 inf>  IrC^o+  iyo  +  C  f  L<*J 

fin/,  ff£.+  iy.) 


(13) 


The  evaluation  of  the  above  expression  for  a  given  cascade  configura¬ 
tion  (i.  e.  ,  a  given  s,  c,  and  at  )  is  fairly  straightforward. 

To  roughly  estimate  the  induced  velocities  in  the  present 
investigation,  we  use  an  even  simpler  approximation  of  representing 
the  blade  by  a  single  vortex  line  of  strength  J~*  ,  located  say  at 
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£  =  0.  For  this  case, 


IT-  L  w 


CotA  f(t.+  i 


the  induced  velocities  v  and  w  are  plotted  against  yD/s  for  various 
values  of  ZQ! ,s  {Fig,  17): 


In  the  present  investigation,  at  the  measuring  radius, 


r  _ 


/.i-2  V* 


where  Va  is  the  mean  axial  velocity.  Also,  c/s,  or  the'ratio  of 
chord  to  spacing,  for  the  rotor  is  1.  27.  If  one  assumes,  for  example, 
that  the  equivalent  single  vortex  is  located  at  the  middle  of  the  blade, 
a  simple  calculation  based  on  the  geometry  of  the  rotor  (Fig.  4)  shows 


that  measuring  Station  6  is  at  a  distance  Zol  s  =  0.  8  from  the 


vortex.  Fig.  17  gives 


<7%r). 


*  '•  *  * 


<£)_„  -  (£) 


I.ol  % 


V* 


0-  7S  2 


*  For  sinusoidal  functions,  the  r.  m.  s.  value  to  the  peak 
value  is  1/\T 2.  The  variation  of  v  and  w  with  yD/s  is,  how¬ 
ever,  not  sinusoidal,  although  not  far  off  from  it. 
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The  measured  value  of  fluctuation  intensity  at  Station  6  is  1.  45%.  Let 
us  assume  that  the  turbulence  intensity  is  the  same  as  its  value  at 
Station  4,  namely  .  30%.  Then  the  fluctuation  intensity  due  to  the 
induced  velocities  alone  would  be 

a/v/U  ^  /■  42.  1 

according  to  Eq.  (5).  This  figure  would  put  the  equivalent  single 
vortex  at  approximately  40%  of  chord  length  from  the  blade  leading 
edge.  This  is  an  entirely  convincing  figure.  Furthermore,  this 
figure  also  checks  the  fluctuation  intensity  at  Station  5.  {A  more  re¬ 
fined  calculation  is  to  use  Eq.  (12)  or  (13).  ) 

If  this  analysis  on  induced  velocities  is  correct  (both  fluctua¬ 
tion  intensity  measurements  and  oscilloscope  photographs  seem  to 
attest  to  its  validity),  then  the  intensity  of  the  nonsteady  induced 
velocities  can  easily  exceed  10%  of  mean  velocity  near  the  rotor 
blade  (see  Fig.  17).  This  would  have  the  tendency  to  accelerate  the 
decay  and  at  the  same  time  increase  the  width  of  the  wakes.  In  the 
recent  investigations  by  Sears  and  Kemp  (Kef.  15)  and  Meyer  (Kef.  7) 
on  the  influence  of  wakes  in  a  turbomachine,  the  wake  width  was  taken 
from  measurements  on  single  airfoils  (Kef.  12)  and  therefore  was 
assumed  to  be  narrow  compared  to  the  chord.  It  appears  that  the 
results  thus  obtained  may  be  somewhat  exaggerated--tbe  actual  wake 
influence,  thanks  to  the  induced  velocities,  is  probably  less  severe. 
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(2)  Trailing  Vortices.  When  the  blade  vortex  strength  varies 
along  the  span,  the  trailing  vortices  shed  from  the  blades  produce  an 
induced  velocity  in  addition  to  that  produced  by  the  bound  vortices  dis¬ 
cussed  previously. 

Using  the  notation  of  Figures  14  and  16,  we  investigate  the 
velocity  at  P  induced  by  the  trailing  vortex  situated  at  x  =  x  and 
along  the  direction  of  f  ,  from  _f  =  f  to  =  oo  ,  Let  us 

express  the  location  of  P  in  the  ^  coordinates  system, 

i.  e.  ,  P  (xq  ,  S  q  ,  ^ 

The  induced  velocity  at  P  due  to  one  trailing  vortex  behind 
one  blade  (the  reference  blade)  is: 


n=-foo  **» 


i 


C2Y  Ol,  +  oL) _ 

(x  -X.)2  +  (rjB+ ns  Cosot. )*■ 


S„  +S  +/}S  °*s 


+  00 

ar 


V  &Z  i  f 

'  /ts-cO'co  -4 


J(  X-Xo^-hCjo  +  nSCosaL)1 

(x-x») 


=====  ■■  Us 

+  (f+s.+nssiocij*) 


dx 

(15) 


<?*  (x-x0f  +(vfo+ns  Ccsnj2- 

I  /  +  _ L-±l±..n*  sin«  - - Idtdx 

l  J(x -*<,)*+ (ye  + ns  coset ) 1  +(f  -#■/,  +  ns  Sin*)*  j  * 


Referring  to  x,  y,  z  coordinates,  the  components  v,  w  of  the 
induced  velocity  are  related  to  q^  in  the  following  way. 

v  =  9  cos  <X 

w  =  9^  Since 

These  velocities  v  and  w  are  to  be  added  to  the  values  given  by 
Eqs.  (9). 

Solutions  to  the  above  equations  have  to  be  worked  out 
numerically  for  each  design. 

It  should  be  mentioned  that  in  analyzing  the  induced  velocities 
due  to  the  rotor,  only  the  influence  of  the  steady  part  of  the  bound 
and  trailing  vortices  of  the  rotor  blades  has  been  considered.  The 
unsteady  part  of  the  bound  and  trailing  vortices  (caused  by  the  inter¬ 
action  of  the  stator)  is  assumed  to  be  of  smaller  order  of  magnitude 
when  compared  to  the  steady  part,  and  hence  can  be  neglected. 
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5d.  Fluctuation  Intensity  Behind  a  Rotor.  To  estimate  the 


intensity  of  the  periodic  fluctuations  caused  by  viscous  wakes  behind 
a  rotor,  it  is  interesting  to  perform  the  following  simple  analysis: 

Let  the  velocity  outside  of  the  wakes  be  unity,  and  the  velocity 
inside  of  the  wakes  be  m.  Therefore,  m  is  a  number  between  0 
and  1,  depending  on  the  distance  downstream  from  the  rotor.  The 
smaller  this  distance,  the  smaller  is  m.  Let  £.  be  the  ratio  of  the 
width  of  the  wake  to  the  spacing  between  two  blades  (Fig.  18).  It  can 
be  easily  shown  that  the  average  velocity  is  1  -  (1  -  m)£  ,  and  the 

r.m.  s.  of  the  fluctuation  is  (1  -  m)  V  £  (1  -  £.  ).  The  intensity  of 
these  periodic  fluctuations  (as  a  ratio  to  the  average  velocity)  is  there¬ 
fore 

-SL  _  -/£('-  u 

U  (-tt)-i  ,161 

The  peak  fluctuation  (in  the  wake  measured  from  the  mean)  is,  again 
as  a  ratio  to  the  average  velocity, 

U"’“*  .  ■  ,1-^  •  (17) 

v  Ct^)-£ 

Fig.  19  shows  Eq.  (16)  plotted  for  the  fluctuation  intensity  in 
terms  of  m  and  £  .  A  close  look  at  the  condition  of  axial  velocity 
fluctuation  at  the  measuring  Station  7  (see  Fig.  10,  Case  7a)  shows 
that,* 


*  Since  a  hot-wire  is  sensitive  only  to  the  velocity  component 
perpendicular  to  it,  in  the  following  u  and  U  will  be  axial 
components  for  Case  7a  and  tangential  components  for  Case  7b. 
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£  =*•  .08 


Uynax/ U  ~  /&  2 


This  gives,  from  Eqs.  (16)  and  (17) 

w  ^  0.8  J^/U  &  5.  5  % 

This  checks  closely  the  measured  value  of  fluctuation  intensity  of 
5.47%,  The  contribution  of  turbulence  to  the  total  fluctuation  in¬ 
tensity  is  probably  small. 

We  now  turn  our  attention  to  Case  7b,  namely  the  tangential 
velocity  fluctuations  behind  the  rotor.  A  simple  analysis  based  on  the 
assumption  that  the  exit  velocity  always  maintains  a  constant  direc¬ 
tion  relative  to  the  blade  (although  the  magnitude  of  the  relative 
velocity  decreases  in  the  wake)  gives  immediately  the  ratio  m'  of 
the  tangential  velocity  in  the  wake  to  the  tangential  velocity  outside 
the  wake  as  a  function  of  m  for  the  axial  velocity  and  the  blade 
geometry.  Referring  to  Fig.  20,  it  can  be  seen  that 

/*  (f~m)  ("w*  -l) 

where,  following  convention,  vu  =  blade  velocity,  =  absolute 
tangential  velocity,  and  Va  =  axial  velocity. 


VADC  TR  55-249 


23 


IX) 


0.  80, 


For  the  case  under  consideration,  = 

then  m'  =  1.  0616.  Eqs.  (16)  and  (17)  can  again  be  applied"  with  m' 
instead  of  m.  It  is  found  that  the  intensity  of  periodic  fluctuation 
due  to  wakes  is 

/jy/u  =  i.  7% 

The  measured  total  fluctuation  intensity  is  2.2%.  From 
Eq.  (5),  the  turbulence  intensity  can  be  estimated, 

Ju}fu  —  J 2. Zz -/.  7*  =  1.4% 


An  examination  of  Fig.  10,  Case  7b,  shows  that  the  turbulence  level 
indeed  appears  to  be  of  the  same  order  of  magnitude  as  the  periodic 
wakes. 

It  might  be  added  that  although  the  accuracy  of  constant- 
current  type  of  hot-wire  anemometer  is  presumably  in  doubt  for 
very  large  fluctuations,  the  value  of  U'max/U  =  18%  obtained  in 
the  wake  of  rotor  blades  checks  very  well  with  the  wake  survey 
behind  inlet  guide  vanes  at  comparable  distance  from  the  blades. 
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43  Blades  Inlet  Guide  Vanes  LGV. 


I-  Cross  Section  of  Testing  Compressor 


Figure  2  Traversing  Apparatus 
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Figure  3  Hot-Wire  Anemometer  Equipment 


Behind  Stator,  Fluctuation 


axial  distance  (inches)  from  trailing  edge  of  inlet  guide  vanes 
Fig.  12.  DECAY  OF  WAKE  BEHIND  INLET  GUIDE  VANES 
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Fig.  13.  FLUCTUATION  INTENSITY  BEHIND  INLET 

GUIDE  VANES 
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